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Almost forty years ago, Carnahan and Starling pub- 
lished in this journal their, now famous, equation of state 
for hard sphere fluid pj. Their derivation was based on 
the simple observation that the leading order virial coef- 
ficients for hard sphere fluid in three dimensions closely 
followed a geometric sequence. The assumption that this 
behavior also extrapolated to higher order virials, allowed 
Carnahan and Starling to explicitly resum the virial ex- 
pansion to find a simple, yet, very accurate equation of 
state. 

Unfortunately, no such accurate equation of state is 
known for the case of lattice gases. This is particu- 
larly frustrating, since lattice models are widely used to 
study many complex fluids ranging from microemulsions 
to electrolytes H S H . In this note, we shall present 
a very simple equation of state which works very well for 
two dimension lattice gas of hard squares and reasonably 
well for three dimension lattice gas of small hard cubes 
at not too high density. 

Our discussion is based on a lattice theory of polymer 
mixtures proposed a long time ago by Flory Q , who de- 
duced the entropy of mixing to be 



S = -fc B [iViln0i + A 2 ln</> 2 ] 



(1) 



where N\ and N 2 are the number of polymers of type 
one and two, while <f>i and (j) 2 are their respective vol- 
ume fractions. The form of Eq.Q is particularly appeal- 
ing since it does not contain any reference to the lattice 
structure and depends only on thermodynamically well 
defined variables. The mixture is assumed to fill all the 
available volume, so that there are no vacancies. If there 
is only one type of polymer occupying a volume fraction 
4>i , the rest of the space is taken to be filled by the solvent 
of 4>2 — 1 — 01 • 

It is clear that the formalism developed by Flory for 
polymer mixtures should be readily applicable to "hard" 
non-attracting lattice gases. Consider, for example, a 
lattice gas of hard hypercubes of volume X d (A is integer 
and the lattice spacing is taken to be 1) on a simple 
hypercubic lattice in d dimensions. The Helmholtz free 
energy of this lattice gas is F m — —TS, since the system 
is athermal. The free energy density is 



/3f m = pln4>+(l-<t>) ln(l 



(2) 
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where p is the particle density and cf> = X d p is the volume 
fraction. 

We note, however, that in the low density limit, Eq. 
(J2J does not reduce to the free energy of the ideal gas 



f3f = p\np- p 



(3) 



Therefore, f m can not be the total free energy of the sys- 
tem, except for the case of A = 1 when Eq. J5J) becomes 
exact. For polymer mixtures, to obtain the total free en- 
ergy, Flory added an extra contribution to Eq.J2J which 
accounted for the conformational degrees of freedom of 
the polymer chains, the so called entropy of disorienta- 
tion 0. This restored the correct low density behavior 
to the theory. For rigid particles, however, the entropy of 
disorientation is identically zero and cannot be the reason 
for the failure of Eq. J2J). 

To recover the correct low density behavior, while 
preserving the simple and thermodynamically appealing 
form of Eq. J5J), we modify f m by introducing a multi- 
plicative factor 17(A) into Eq. (J2J, 



/3f = p]n4> + g(\)(l-<l>) ln(l 



(4) 



This equation can be interpreted as an interpolation be- 
tween the low density limit governed by the particles, 
and the high density limit in which defects, the "holes" , 
become relevant. The total number of holes, however, is 
not fixed since the vacancies can change their size and 
shape, so that the prefactor appearing in front of the 
second term of Eq. |0J is the effective number density. 

The requirement that in the low density limit Eq. Q 
must reduce to Eq. J3J), uniquely determines the func- 
tional form of <?(A) yielding 



<7(A) 



1 + dlnA 
A^ 



(5) 



Note that (7(1) = 1, so that Eq. |0J reduces to the exact 
free energy for a lattice gas of A = 1. The chemical po- 
tential within the modified Flory approximation (MFA) 



is 



0H = -(1 + din A) [ln(l -<P) + 1) +In0+1 



(G) 



In Fig. n we compare the value of the chemical po- 
tential obtained within the MFA with the results of the 
Monte Carlo simulations for a gas of hard squares of dif- 
ferent sizes A. The simulations were performed using 
the grand-canonical ensemble at fixed volume V, tem- 
perature T, and chemical potential /1, with trial moves 
insertion and removal of particles as well as attempts to 
diffuse 8]. The agreement is excellent for all As tested, 
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FIG. 1: Chemical potential versus volume fraction for various 
two dimensional lattice gases. The symbols are the MC results 
while the lines are the predictions of the MFA Eq. @, with 
d = 2. 
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FIG. 2: Chemical potential versus volume fraction for A = 2 
lattice gas. The symbols are the MC results; the solid curve is 
the equation of state obtained in this Note, Eq. (|HJ ; the dashed 
line is the equation of state obtained using the Fundamental 
Measure Density Functional Theory, Ref. |Tll| . which is also 
the same as the one found earlier by Tcmperley 10]. The 
order-disorder transition occurs at <j) « 0.93. 



In Figs. |3 and 0] we also show the equations of state 
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FIG. 3: Chemical potential versus volume fraction for a gas 
of hard rods with A = 1, 2 and 3 in d = 1. The symbols are 
the exact value of the chemical potential [T^| while the lines 
are the predictions of the MFA Eq. © with d = 1. 
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FIG. 4: Chemical potential versus volume fraction for various 
three dimensional hard cube lattice gases. The points are the 
MC results while the lines are the predictions of the MFA Eq. 
© with d = 3. 



up to quite high volume fractions. However, similarly to 
the Carnahan-Starling equation of state, MFA also fails 
to notice the phase transition between the disordered and 
ordered (columnar) phases present at high volume frac- 
tions In Fig. |21 we compare the accuracy of the MFA 
with the earlier equation of state derived by Temper- 
ley 01 which is also identical to the one recently found 
using the Fundamental Measure theory [Tl| . 

It is curious that all the MC curves for different val- 
ues of A intersect at approximately the same point. This 
property is also captured by the MFA, which predicts 
that all the chemical potentials for different As are equal 
when the volume fraction satisfies ln(l — X ) = —1, in- 
dependent of d. The value (j> x = 0.632121 is in excellent 
agreement with the intersection point observed in the 
Monte Carlo simulations. 



for d = 1 hard rods, and d = 3 hard cubes. In the case of 
d = 1 the exact free energy is known |l2( . Although still 
quite good, the agreement between the simulations and 
the MFA deteriorates more rapidly with increasing A for 
d = 1 and d = 3 than for d = 2. 

The high degree of accuracy of the MFA in d — 2 is 
quite surprising in view of the crudeness of the approx- 
imation. It also suggests that there should be a more 
direct way to arrive at Eq. or some other such equa- 
tion of state 0, , generally valid for non-attracting 
lattice gases of arbitrary A. In the absence of such general 
theory, the very simple Eq. Q should be useful for con- 
structing lattice mean-field theories for various complex 
systems. 

This work was supported in part by the Brazilian sci- 
ence agencies CNPq, CAPES, and FAPERGS. 
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